Spin pumping and magnetization dynamics in metallic multilayers 
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We study the magnetization dynamics in thin ferromagnetic films and small ferromagnetic parti- 
cles in contact with paramagnetic conductors. A moving magnetization vector causes "pumping" of 
spins into adjacent nonmagnetic layers. This spin transfer affects the magnetization dynamics simi- 
lar to the Landau-Lifshitz-Gilbert phenomenology. The additional Gilbert damping is significant for 
small ferromagnets, when the nonmagnetic layers efficiently relax the injected spins, but the effect 
is reduced when a spin accumulation build-up in the normal metal opposes the spin pumping. The 
damping enhancement is governed by (and, in turn, can be used to measure) the mixing conductance 
or spin-torque parameter of the ferromagnet-normal-metal interface. Our theoretical findings are 
confirmed by agreement with recent experiments in a variety of multilayer systems. 



PACS numbers: 72.25.Mk,75.70.Cn,76.50.-|-g,76.60.Es 



I. INTRODUCTION 

Spin-polarized transport through magnetic multilay- 
ers is the physical origin of many exciting phenom- 
ena such as giant magnetoresiatajice and spin-current- 
induced magnetization reversal.aal3 It has attracted at- 
tention in the basic physics community and industry over 
the last decades, but there are still open fundamental 
questions. So far, the main research activity has been 
focused on the dc transport properties of these systems. 

Ac magnetotransport has attracted considerably less 
attention than its dc counterpart. In a recent paper,cl 
we reported a novel mechanism by which a precessing 
ferromagnet "pumps" a spin current into adjacent non- 
magnetic conductors proportional to the precession fre- 
quency, using a formalism analogous to that for the,adi- 
abatic pumping of charges in mesoscopic systems.El We 
showed that spin pumping profoundly affects the dynam- 
ics of nanoscale ferromagnets and thin films, by renormal- 
izing fundamental parameters such as the gyromagnetic 
ratio and Gilbert damping parameter, in agreement with 
experiments.^ 

The switching characteristics of a magnetic system 
depends in an essential way on the Gilbert damping 
constant a. In magnetic field-induced switching pro- 
cesses, for example, a governs the technologically impor- 
tant magnetization reversal tipjie of ferromagnetic parti- 
cles. Its typical intrinsic valucQ < 10~^ for transition 
metal ferromagnets is smaller than its optimal value of 
a > 10~^ for the fastest switching rates. Q The present 
mechanism allows engineering of the damping constant 
by adding passive nonmagnetic components to the sys- 
tem and/or adjusting the geometry to control spin flow 
and relaxation rates described in this paper, thus helping 
to create high-speed magnetoelectronic devices. Also, in 
spin-current-induced magnetization reversal, the critical 
switching current is proportional to aJ3 

For some time it has been understood that a 



ferromagnet-normal-metal {F-N) interface leads to a dy- 
namical coupling between the ferromagnetic magnetiza- 
tion and the spin^ xiJie. conduction-band electrons in 
the normal metalE'ElElEjO'tJ More recently, several the- 
oretical frameworks were put forward proposing a mech- 
anism for .maguetization damping due to F-N interfacial 



processes.clEirEil This F-N coupling becomes important 
in the limit of ultrathin (<10 nm) ferromagnetic films 
and can lead to a sizable enhancement of the damping 
constant. 

Our theory is based on a new physical picture, accord- 
ing to which the ferromagnetic damping can be under- 
stood as an adiabatic pumping of spins into the adjacent 
normal metals. u This spin transfer is governed by the re- 
flection and transmission matrices of the system, analo- 
gous to the scattering theory of transport and interlayer 
exchange coupling. The microscopic expression for the 
enhanced Gilbert damping and the renormalized gyro- 
magnetic ratio can be calculated by simple models or by 
flrst-principles band-structure calculations without ad- 
justable parameters. The present theory therefore allows 
quantitative predictions for the magnetization damping 
in hybrid systems that can be tested by experiments. 

The Gilbert damping constant in Jjl iiTi fprjon iagnetic 
fllms has been experimentally studicdS'liirlijOO by mea- 
suring ferromagnetic-resonance (FMR) linewidths. In 
the regime of ultrathin ferromagnetic films, a was in 
some cases found to be quite large in comparison with 
the bulk value ao, and sensitively depending on the sub- 
strate and capping layer materials. For example, when a 
20-A-thick permalloy (Py) film was sandwiched between 
two Pt layers, its damping was found to be a ~ 10~^, 
but recoitered its bulk value a ~ 10^^ with a Cu buffer 
and cap.l3 Heinrich et alS3 observed an enhanced damp- 
ing of <20 A-thick Fe films when they were grown on Ag 
bulk substrates but no significant change in the damping 
constant was seen for films grown on GaAs even when 
the film thickness was reduced down to several atomic 
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monolayers. t2l We will demonstrate here that our theory 
explains all these experimental findings well. 

Previously, we studied the situation when the normal- 
metal layers adjacent to the ferromagnetic films are 
perfect spin sinks, so thatrnthe spin accumulation in 
the normal metal vanishes.Q Here this theory is gen- 
eralized to consider the spin accumulation which en- 
ables uSpijO .explain experimental findings for various F-N 
systems0oll3 in a unified framework based on the spin- 
pumping picture. 

The paper is organized as follow. In Sec. |l|and the Ap- 
pendix, the basic formalism of the adiabatic spin-pump 
theoryQ is derived using a scattering-matrix approach, 
and an alternative derivation is given for finite systems, 
which is based on the con servation of energy and angular 
momentum. In Sees. Ill and IV, we solve the diffusion 



equation to describe transport of injected spins in single 
and composite normal-metal layers. The spin loss due 
to spin-orbit or other spin- flip processes is accounted for, 
leading to an overall damping of the ferromagnetic mag- 
netization precession. In particular, we use the theory 
to analyze the representative case of Gilbert damping in 
Py-Pt, Py-Cu, and Py-Cu-Pt hybrids, showing an excel- 
lent agreement between our theory and the experimental 
results.l3ll3 The last Sec. ^ is devoted to discussions and 
conclusions. 



II. PRECESSION-INDUCED SPIN PUMPING 

Consider an N-F-N junction as in Fig. |l| and Fig. | m 
the Appendix. Without a voltage bias, no spin or charge 
currents flow when the magnetization of the ferromagnet 
is constant in time. When the magnetization direction 
starts precessing (as, e.g., under the influence of an ap- 
plied magnetic field), a spin current IP"™p is pumped out 
of the ferromagnet This current into a given N layer 
depends on the complex- valued parameter A = Ar + iAi 
(the "spin-pumping conductance" ) by 
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AjVci X — Ai—— 

dt dt 



(1) 



Here, the time-dependent order parameter of the ferro- 
magnet is a unit vector m(t), assuming a monodomain 
magnet with a spatially uniform magnetization at all 
times. A detailed derivation of Eq. (|^) based on the 
scattering-matrix theory of transport is presented in the 
Appendix. A — g^^ — depends on the scattering ma- 
trix of the ferromagnetic film since 

g'^'^'^^[5„„-C„(C„)*] (2) 

mn 

is the dimensionless dc conductance matrix0lll and 



Here (see Fig. ||), rj^„ (fmn) is a reflection coefficient for 
spin-up (spin-down) electrons on the normal-metal side 
and tj]j„ {t'mn) is a transmission coefficient for spin-up 
(spin-down) electrons across the ferromagnetic film from 
the opposite reservoir into the normal-metal layer, where 
m and n label the transverse modes at the Fermi energy 
in the normal-metal films. Note that the magnetization 
can take arbitrary directions; in particular, m(t) may be 
far away from its equilibrium value. In such a case, the 
scattering matrix itself can depend on the orientation of 
the magnetization, and one has to use A{m) in Eq. (jl|). 

When the ferromagnetic film is thicker than its trans- 
verse spin- coherence length, d > 7r/(fcp — kp), where 

/cp^^-* areulhe spin-dependent Fermi wave vectors, t^^ 
vanishes,L2l the spin pumping through a given F-N in- 
terface is governed entirely by the interfacial mixing con- 
ductance A — q^^ = n'^^ 
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n 



and we can consider 



only one of the two interfaces. This is the regime we are 
focusing on in this paper. Note that the conductance 
matrix 5°''^ defined in Eq. (|^) has to be renormalized 
for highly transparent interfaces in columnar geometries 
(by properly subtracting Sharvin resistance contributions 
from the inverse conductance parameters), as discussed 
in Ref. 
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FIG. 1: A ferromagnetic film F sandwiclied between two non- 
magnetic reservoirs A'^. For simplicity of tfie discussion in tliis 
section, we mainly focus on tfie dynamics in one (rigfit) reser- 
voir while suppressing the other (left), e.g., assuming it is 
insulating. The spin-pumping current IP"™p g^j^^j ^j^g gpjjj g^^_ 
cumulation fi^ in the right reservoir can be found by conser- 
vation of energy, angular momentum, and by applying circuit 
theory to the steady state ir""" = Is''*- 

As shown before,^ the spin current [Eq. (|l|)] leads to 
a damping of the ferromagnetic precession, resulting in 
a faster alignment of the magnetization with the (effec- 
tive) applied magnetic field Hoff. In the derivation by 
the time-dependent scattering theory, the pumped spins 
are entirely absorbed by the attached ideal reservoirs. In 
the following, it is shown that Eq. (pi) can be also de- 
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rived for a finite system by observing that the enhanced 
rate of damping is accompanied by an energy flow out 
of the ferromagnet, until a steady state is established 
in the combined F-N system. For simplicity, assume a 
magnetization which at time t starts rotating around the 
vector of the magnetic field, m{t) _L Hcff. In a short in- 
terval of time St, it slowly (i.e., adiabatically) changes to 
m{t + St) — m{t) + Sm. In the presence of a large but fi- 
nite nonmagnetic reservoir without any spin-flip scatter- 
ing attached to one side of the ferromagnet, this process 
can be expected to induce a (small) nonvanishing spin 
accumulation 



/X, = / deTi[af{e)] 



(4) 



where & is the Pauli matrix vector and /(e) is the 2x2 
matrix distribution function at a given energy e in the 
reservoir .M For a slow enough variation of m{t), this 
nonequilibrium spin imbalance must flow back into the 
ferromagnet, canceling any spin current generated by the 
magnetization rotation, since, due to the adiabatic as- 
sumption, the system is always in a steady state. 

Let us assume for the moment that the spins are 
accumulated in the reservoir along the magnetic field, 
fig II Heff. Flow of Ns spins into the normal metal 
transfers energy AEn — NsfJ-s/'^ and angular-momentum 
ALn = Nsfi/2 (directed along Hcff). By the conser- 
vation laws, AEp = —AEm and ALp = —AL^, for 
the corresponding values in the ferromagnet. Using the 
magnetic energy AEp = jALpHcs, where 7 is the ab- 
solute gyromagnetic ratio of the ferromagnet, we then 
find that NsHsf^ = jNs{h/2)HcS- It then follows that 
fis = hrfHcS — huj, where w — jHcS is the Larmor fre- 
quency of precession in the effective field: The spin-up 
and spin-down chemical potentials in the normal metal 
are split by /i^ = /zw, the energy corresponding to the 
frequency of the perturbation. For a finite angle 9 be- 
tween fig and Hcff, the same reasoning would lead to 
fj,s = fko COS 9, which is smaller than the "energy boost" 
huj of the time-dependent perturbation, thus justifying 
our initial guess. l— . 

We can employ magnetoelectronic circuit theorylHl to 
derive an expression for the backflow of spin current 1^'*'='^ 
which, as argued above, has to be equal to the pumping 
current 1?""^? = I^a'^'^: 
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dm 

~dt 



(5) 



Here, we used fjLg ^ fiio and /i.^ _L m, since by the conser- 
vation of angular momentum, the spin transfer is propor- 
tional to the change in the direction Sm _L m. We thus 
recover Eq. (|^) for the case of a single and finite reser- 
voir. It is easy to repeat the proof for an arbitrary initial 
alignment of m{t) with Hcff. Furthermore, a straightfor- 
ward generalization of this discussion to the case of the 



N-F-N sandwich structure recovers our previous result 
[Eq. (0)]. 

The expressions for the adiabatic spin pumping are not 
the whole story, since spin-flip scattering is an important 
fact of life in magnetoelectronics. In Ref. ^, we only 
considered the extreme situation where the normal-metal 
layer is a perfect spin sink, so that all spins injected by 
jpump ],gjg^-5j- j^y spin-flip processes or leave the system; 
the total spin current through the contact was, therefore, 
approximated by Is « ip"™p and I^^'^'^ « 0. Here, we 
generalize that treatment to self-consistently take into 
account the spin build-up in the normal metal at dynamic 
equilibrium. We then find the contribution to due to 
the spin-accumulation-driven current I^^'^'^ back into the 
ferromagnet: 



I — IP" 

J-.S J-c 



rback 



(6) 



which vanishes in the absence of spin-flip scattering. 

The spin current out of the ferromagnet carries angular 
momentum perpendicular to the magnetization direction. 
By conservation of angular momentum, the spins ejected 
by Is correspond to a torque r = — I^ on the ferromagnet. 
If possible interfacial spin-flip processes are disregarded, 
the torque t is entirely transferred to the coherent mag- 
netization precession. The dynamics of the ferromagnet 
can then be described b¥r-a| generalized Landau-Lifshitz- 
Gilbert (LLG) equatioiJI-El 

dm dm 7 ^ 

— = -7m X Hcff + aom X — + ^1. , (7) 

where ao is the dimensionless bulk Gilbert damping con- 
stant, Mg is the saturation magnetization of the ferro- 
magnet, and V is its volume. The intrinsic bulk constant 
q;o is smaller than the total Gilbert damping a = ao + Q^'- 
The additional damping a' caused by the spin pumping 
is observable in, for example, FMR spectra and is the 
main object of interest here. 



III. SPIN-ACCUMULATION-DRIVEN 
BACKFLOW IN THE F-N AND N-F-N 
MULTILAYERS 

The precession of the magnetization does not cause any 
charge current in the system. The spin accumulation or 
nonequilibrium chemical potential imbalance fjbg {x) [sim- 
ilar to Eq. (Q), but spatially dependent now] in the nor- 
mal metal is a vector, which depends on the distance 
from the interface < a; < L, where L is the thickness 
of the normal-metal film, see Fig. ^. 

When the ferromagnetic magnetization steadily ro- 
tates around the z axis, m x m and the normal-metal 
spin accumulation Hs{x) are oriented along z, as depicted 
m Fig. |. There is no spin imbalance in the ferromag- 
net, because fig is perpendicular to the magnetization 
direction m. As shown below, the time-dependent fig is 
also perpendicular to m even in the case of a precessing 
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FIG. 2; Schematic view of the F-N bilayer. Precession of the 
magnetization direction m.(t) of the ferromagnet F pumps 
spins into the adjacent normal-metal layer A'' by inducing a 
spin current 1^"™^. This leads to a build-up of the normal- 
metal spin accumulation which either relaxes by spin-flip scat- 
tering or flows back into the ferromagnet as I^^'^'^. In contrast 
to Fig. ^, the N layer here is not an ideal reservoir but rather 
a film of the same cross section as the magnetic layer F\ the 
spin accumulation is position [x) dependent. 



Lo SO that K « . For a static applied field 

of 1 T, typically oj ^ 10^^ s^^. The elastic scattering 
rate corresponding to a mean free path of Ad ^ 10 nm 
is T~[^ ~ 10^^ s~^. Consequently, the derivation be- 
low is restricted to metals with a ratio of spin-conserved 
to spin-flip scattering times e = Td/rsf > 10^'^. In 
practice]^ this condition is easily satisfied|_with higher 
impurity atomic numbers Z (as e scales a.aEj Z*). The 
high-frequency limit uj > '''7/ 1 on the other hand, is rel- 
evant for hybrids with little spin-flip scattering in the 
normal metal, and was discussed in the context of the 
spin-battery concept .c3 Nevertheless, we will see that a 
sizable Gilbert damping enhancement requires a large 
spin-flip probability e > 10^^ (thereby guaranteeing that 
u) <C T^f') unless the frequency is comparable with the 
elastic scattering rate in the normal metal. The latter 
regime will not be treated in this paper. 

Using relation D — VpTc\/3 between the diffusion co- 
efficient D (in three dimensions), the Fermi velocity up, 
and the elastic scattering time Td, we find for the spin- 
diffusion length 



Asd = "^F V TclTsf/3 . 



(11) 



ferromagnet with time-dependent instantaneous rotation 
axis, as long as the precession frequency uj is smaller than 
the spin-flip rate t^^ in the normal metal. 

The spin accumulation diffuses into the normal metal 

as 



(8) 



where D is the diffusion coefficient. The boundary condi- 
tions are determined by the continuity of the spin current 
from the ferromagnet into the normal metal at x = and 
the vanishing of the spin current at the outer boundary 
X — L: 



X 
X 



: d.^n, 
L : dxUs 



'2{hMSD)-% , 



(9) 



where Af is the (one-spin) density of states in the film and 
S is the area of the interface. The solution to Eq. (||) with 
the boundary conditions [Eqs. (||)] is 



cosh k{x — L) 2Is 
sinh kL TtMSDhl 



(10) 



with the wave vector k — A^ vT^flaJr^, where Asd = 
^/ DTsi is the s pin -flip diffusion length in the normal 
metal. In Ref. E2 we used arguments similar to those 
in the present paper to calculate the spin accumulation 
( p^ ) generated by the precessing magnetization. While 
the size of the effect and its relevance for spintronic ap- 
plications are detailed in Ref. |2^, in this work we focus 
on the role of the spin accumulation in the dynamics of 
the ferromagnetic magnetization. 

We assume in the following that the precession fre- 
quency Lo is smaller than the spin-flip relaxation rate 



An effective energy-level spacing of the states participat- 
ing in the spin-flip scattering events in a thick film can 
be defined by 



{AfSK 



(12) 



The spin-accumulation-driven spin current 1^'^'=^ through 
the interface readstJ 



fback 



1 

8^ 



2gTi;x,(x = 0) + 2gl'ra x ^l^{x = 0) 



-2.9Ti)(m.M,(x = 0))m] . (13) 



Substituting Eq. ( |Io|) into Eq. (|T^, we find for the total 
spin current [Eq. fflj] 



I IP 

J-.s — J-e 



5 



TT 
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2gl^ls + 2gl^m x I, 
-2gl^) (m-I,)m] , 



(14) 



where the spin current returning into the ferromagnet is 
governed by the "backflow" factor /3, 



TsfSsd/h 



tanh(_L/A; 



sd) 



(15) 



When the normal metal is shorter than the spin-diffusion 
length (i < Asd), l3 Ts[S/h, where S = (7V5L)"^ 
is the energy-level splitting. In the opposite regime of 
thick normal metals {L 3> Asd), P Tgidsd/h. Basically, 
/3 [Eq. (|l^)] is therefore the ratio between the energy 
level spacing of the normal-metal film with a thickness 
Lsf — min(L, Asd) and the spin-flip rate. 

By inverting Eq. (|lj), we may express the total spin 
current in terms of the pumped spin current I?"™? 
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[Eq. (0)] 



1 + figl^ + 



n\2 



1 + pgl^ 



1 + /3.g, 



n 



mx I . 



After substituting Eq. (g) into Eq. (|l|), we recover the 
form of Eq. for the total spin current , but with a 
redefined spin-pumping conductance A = Ar+ iAi 



h f ~ dm ~ dm 

1= = — Aj in X Ai 

47r V dt dt 



A can be expressed in terms of the mixing conductance 
g^'^ and the backflow factor /? by 



A,, 



1 pglHi + PgV)-^ 

-Pgl^{l + (3gl^)-^ 1 



1 + 



Ti 



M\2 



\ + (3g 



Ti 



9. 
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angular momentum that diffuses back into the ferromag- 
net. It was disregarded in Ref. ^ where the normal metal 
was viewed as a perfect spin sink. Because spins accumu- 
late in the normal metal perpendicular to the ferromag- 
netic magnetization, the spin-accumulation-driven trans- 
port across the F-N contact, as well as the spin pumping, 
(16) is governed by a mixing conductance. This explains why 
the other components of the conductance matrix [Eq. (H)] 
do not enter Eq. (|l9|) . 

We now estimate the numerical values of the param- 
eters in Eq. ( p^ ) for transition metal ferromagnets Fe, 
Co, and Ni, in contact with relatively clean simple nor- 
mal metals Al, Cr, Cu, Pd, Ag, Ta, Pt, and Au. For an 
isotropic electron gas, M — /cp/(7r/iWF). Using Eqs. ( [Tl| ) 
and (p^, we find h/{5sdTsi) = 4:^/e/3Nc-h, where Nch = 
Skp I (47r) is the number of transverse channels in the nor- 
mal metal and e = Td/Tsfis the spin- flip probability at 
each scattering. In Ref. g6[ g^^ was calculated for Co- 
Cu and Fe-Cr interfaces by first-principles band-structure 
calculations. It was found that irrespective of the interfa- 
cial disorder, g^J- w TVch for these material combinations. 
As shown in Ref. |l^ has to be renormalized in such 
(18) limit, making the effective conductances about twice as 
large. We thus arrive at an estimate 



(17) 



It has been showncS that for realistic F-N interfaces 

(The latter approximation 



g]'^ < gl'', so that g 
will be implied for the rest of the paper.) In this impor- 
tant regime, Ai vanishes and the term proportional to Ar 
m Eq. Int) has the same form as and therefore enhances 
the phcnomenological Gilbert damping. This can be eas- 
ily seen after substituting Eq. ( [TtI ) into Eq. (|^) : The last 
term on the right-hand side of Eq. (0) can be combined 
with the second term by defining the total Gilbert damp- 
ing coefficient a = ao + a' , where 



Ti 



'9: 



Ti 



1 + 5 



Ti. 



TsfSsd/h 



tanh(L/Asd) 



-1 Ti 
9L9 ' ^ 



Anil 



(19) 



is the additional damping constant due to the interfacial 
F-N coupling. Here, gL is the g factor and ji is the total 
film magnetic moment in units of /iB- Equation. ([19) is 
the main result of this section. When L ^ oo, Eq. 
reduces to a simple result: a' = gLglg/ {4:Tth), where 



1 _ 1 



911, 9 



Ti 



(20) 



Here Rgd = TsfSsd/h is the resistance (per spin, in units of 
/i/e^) of the normal-metal layer of thickness Xgd- [Which 
follows from the Einstein's relation a = e^DAf connect- 
ing conductivity a with the diffusion coefficient D, and 
using Eq. (^.] It follows that the effective spin pump- 
ing out of the ferromagnet is governed by gl^, i.e., the 
conductance of the F-N interface in series with diffusive 
normal- metal film with thickness Asd-lij 

The prefactor on the right-hand sideuof Eq. ( [l9| ) sup- 
presses the additional Gilbert dampingcl due to the spin 



1+ [\/itanh(L/A,d)] ' 



(21) 



where a'^ = g^g^^ / [Att^) is the Gilbert damping en- 
hancement assuming infinite spin-flip rate in the noepial 
metal Tst 0, i.e., treating it as a perfect spin sinkfl 

It follows that only for a high spin-flip probability 
e > 10^^, the normal- metal film can be a good spin sink 



so that 



a'^. This makes the lighter metals, such 



as Al, Cr, and Cu, as well as heavier metals with only s 
electrons in the conduction band, such as Ag, Au, and 
Ta less effective spin sinks since these metals have a rela- 
tively small spi fl -odpit coupling, typically corresponding 
to e < 10-2 .EEZlEa Heavier elements with Z > 50 and 
p or d electrons in the conduction band, such as Pd and 
Pt, on the other hand, can be good or nearly perfect 
spin sinks as they have a much larger e > 10-^0 This 
conclusion explains the hierarchy of the observed Gilbert 
damping enhancement in Ref. ^: Pt has about 2 electrons 
per atom in the conduction band, which are hybridized 
with d orbitals, and a large atomic number Z — 7% and, 
consequently, leads to a large magnetization damping en- 
hancement in the N-F-N sandwich for thin ferromagnetic 
films. Pd which is above Pt in the periodic table having 
similar atomic configuration but smaller atomic number 
Z — A6 leads to a sizable damping, but smaller than for 
Pt by a factor of 2. Ta is a heavy element, Z = 73, but 
has only s electrons and the damping enhancement is an 
order of magnitude smaller than in Pt. Finally, Cu is a 
relatively light element, Z = 29, with s electrons only 
and does not cause an observable damping enhancement 
at all. According to Eq. (^ij), a sufficiently thick active 
layer, L > Xsd, is also required for a sizable spin relax- 
ation. 
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The limit of a large ratio of spin-flip to non-spin- 
flip scattering e ^ 1 deserves special attention. In this 
regime, Eq. (|2^) does not hold, since by using the diffu- 
sion equation [Eq. (||)] and boundary conditions [Eqs. (|^)] 
we implicitly assumed that e ^ 1. If e > lO"'^, on 
the other hand, even interfacial scattering alone can ef- 
ficiently relax the spin imbalance, and such films, there- 
fore, are good or nearly perfect spin sinks (so that 
a' Q^'oo)i regardless of their thickness (in particular, 
they can be thinner than the elastic mean free path). 

Infinite vs vanishing spin-fiip rates in the normal metal 
are two extreme regimes for the magnetization dynam- 
ics in F-N bilayers. In the former case, the damping 
constant a = ao + gLg^^ /{4:Ttij,) is significantly enhanced 
for thin ferromagnetic films, whereas in the latter case, 
a = ao is independent of the ferromagnetic film thick- 
ness. Experimentally, the two regimes are accessible by 
using Pt as a perfect or Cu as a poor spin sink in con- 
tact with a ferromagnetic thin film, as done in Ref. ^ for 
N-Py-N sandwiches. (Using the N-F-N trilayer simply 
increases a' by the factor of 2, as compared to the F-N 
bilayer, due to the spin pumping through the two inter- 
faces.) The measured damping parameter G = ^M^a is 
shown in Fig. |^ by circles. 
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FIG. 3: Circles show measured (Ref. |6|) Gilbert parameter G 
of a permalloy film with thickness d sandwiched between two 
normal-metal (Pt or Cu) layers. Solid lines are predictions 
of our theory with two fitting parameters. Go and g^^-Py 
bulk damping and Py-Pt mixing conductance, respectively, 
see Eq. (p2). 

For the Cu-Py-Cu trilayer, our theory predicts G{d) — 
Go, while for the Pt-Py-Pt sandwich 



G{d) = Go 



{gL^^Brg^^S-' 



as a function of ferromagnetic film thickness d. The 
Py g factor is gL ~ 2.1B These expression agree with 
the experiments for Go = 1-0 x 10® s^^ and g^^S^^ = 
2.6 X 10^^ cm~^ (see Fig. ||). Both numbers are very rea- 
sonable: Go equals the bulk value 0.7— 1.0 x 10® s~^ for 
Py,c3 while g^^S^^ compares well with g^^S^^ « 1.6 x 
10^^ cm~^ found in angular-magnetoresistance (aMR) 
measurements in Py-Cu hybrids.Ej (We recall that here 
one has to use the renormalized mixing conductance g^^, 
in the notation of Ref. |l9[) In fact, since Pt has two con- 
duction electrons per atom, while Cu-only one, and they 
have similar crystal structures, we expect g^^ to be larger 
in the case of the Py-Pt hybrid, justifying the value used 
to fit the experimental data. We have thus demonstrated 
that the additional damping in ferromagnetic thin films 
can be used to measure the mixing conductance of the 
F-N interface. 



IV. MAGNETIC DAMPING IN F-N1-N2 
TRILAYER 



In this section we consider ferromagnetic spin pumping 
into a bilayer N1-N2 normal-metal system, see Fig. |[ It 
is assumed that the spins are driven into the first normal- 
metal film {Nl) of thickness L. While in Nl, spins are 
allowed to diffuse through the film, where they can relax, 
diffuse back into the ferromagnet, or reach the second 
normal-metal layer {N2). N2 is taken to be a perfect 
spin sink: spins reaching N2 either relax immediately 
by spin-fiip processes or are carried away before diffusing 
back into Nl. We show that measuring the ferromag- 
netic magnetization damping as a function of L in this 
configuration can be used to study the dc mixing con- 
ductance of the two Nl film interfaces as well as the Nl 
spin-diffusion time. 
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FIG. 4: Same as Fig. ^, but now the normal-metal system 
is composed of a bilayer N1-N2. Ferromagnetic precession 
pumps spins into the first normal-metal layer Nl. The spin 
build-up in Nl may flow back into the ferromagnet F as spin 



current I, 



relax in Nl, or return to the second normal- 



27rft 



(22) 



metal layer N2 as spin current Igf^ . The spin accumulation 
in N2 is disregarded since the layer is assumed to be a perfect 
spin sink. 
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The analysis in this sectimi was inspired by exper- 
iments of Mizukami et aZ.,t3 who in a follow-up to 
their systematic study of Gilbert damping in N-Py-N 
sandwicheSjH studied magnetization damping in Py-Cu 
and Py-Cu-Pt hybrids as a function of Cu film thickness 
L. The measured damping parameter G is shown by cir- 
cles in Fig. ^. As shown in the preceding section, Cu is 
a poor sink for the pumped spins, while Pt is nearly a 
perfect spin absorber, thus identifying the Cu film with 
Nl and the Pt layer with N2. 

We use the same notation as in the previous section 
to discuss the F-Nl spin pumping with subsequent spin 
diffusion through Nl. Similar to Eqs. (M), the boundary 
conditions for the diffusion equation (|^) in the normal 
metal Nl are now: 



(23) 



Isi and ls2 are the total spin currents through the left 
(x — 0) and right {x — L) interfaces, respectively. I^i 
(similarly to [Eq. (^)] in the previous section) in- 
cludes the pumped spin current [Eq. (|l|)] and the spin- 
accumulation-driven spin current [Eq. (p^)] contribu- 
tions. Is2, on the other hand, is entirely governed by 
the N1^N2 spin-accumulation-driven flow 



(24) 



where g is the conductance per spin of the N1-N2 inter- 
face. 

Solving the diffusion equation (||) with the boundary 
conditions (23), we find the spin current I^i as we did 
in the preceding section. The Gilbert damping enhance- 
ment due to the spin relaxation in the composite normal- 
metal system is then given by 



1 + 



1^ Tsf(5sd 1 + tanh(£/Asd)gTsf^sd/fe 
h tanh(L/Asd) + .gTsf(5sd//i 



-1 Ti 



Anil 
(25) 

Setting 5 = decouples the two normal-metal systems 
and reduces Eq. (|5|) to Eq. ( |l9|) giving the damping co- 
efficient of the F-Nl bilayer. From Eq. (^5|), we get for 
the Py-Cu(L)-Pt trilayer 



G{L) = Go + 



TsfSsd 1 + tanh(L/Asd)gTsAd/fe 



tanh(L/Asd) + gTsfSsd/h 



2h d 
and for the Py-Cu(i) bilayer (putting 

-1 



5 = 0) 



G{L) = Go 



1 



tanh(L/Asd) 



2h d 



(27) 

In the experiments, the permalloy thickness d = 30 A is 
fixed and the Cu film thickness L is varied between 3 and 



1500 nm as show n by the circles in Fig. ||. Our theoret- 
ical results, Eqs. ( pq) and (p7|), are plotted in Fig. ^ by 
solid lines. We use the following parameters: The bulk 
dampingE^ Go = 0.7 x 10® s^^; the spin-flip probability 
e — 1/700 and the spin-diffusion length Asd = 250 nm 
for Cu (which correspond to elastic mean free path Ad = 
VSeAsd — Ifi |PJ))| in agreement with values reported 
in literaturepSEj fti^S'^i ^ 1.6 x 10^^ cm~2 from the 
aMR measurements Jl3 and gS^^ — 3.5 x 10^^ cm~^ for 
the Cu-Pt contact, which lies between values for the ma- 
jority and minority carriers as measured and calculatedEil 
for the Cu-Co interface. Figure || shows a remarkable 
agreement (within the experimental error) between the 
measurements and our theory. It is important to stress 
that while the profiles of the trends displayed in Fig. ^ re- 
veal the diffusive nature of spin transfer in the Cu spacer, 
they cannot be used to judge the validity of a detailed 
mechanism for spin injection (relaxation) at the Py-Cu 
(Cu-Pt) interface. The case of our spin pumping picture 
is strongly supported by the normalization of the curves 
(in agreement with experiment), which are governed in 
our theory by quantities known from other sources. 



# Py-Cu(L)-Pt 
O Py-Cu(L) 




(26) (§) 



FIG. 5: Circles show the measurements by Mizukami et al. 
(Ref. ^) of the Gilbert damping in Py-Cu-Pt trilayer and Py- 
Cu bilayer as a function of the Cu buffer thickness L. Solid 
lines are our theoretical prediction according to Eqs. (Eq) and 



The trends in Fig. ^ can be understood as follows. 
Since Cu is a poor spin sink, a Py-Cu contact with a 
single Cu film does not lead to a significant damping 
enhancement. The small spin-fiip ratio, e 1, causes 
most of the spins transferred into the normal-metal layer 
to be scattered back and relax in the ferromagnet be- 
fore flipping their direction in the Cu buffer. This leads 
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only to a small damping enhancement, which saturates 
at L 3> Asd and vanishes in the limit L <C Agd- The 
situation changes after a Pt film, a very good spin sink, 
is connected to the bilayer: If the normal-metal layer is 
smaller than the elastic mean free path, L <^ Ad, the spin 
accumulation is uniform throughout the Cu buffer. The 
spin pumping will now be partitioned. A fraction of the 
pumped spins reflects back into the ferromagnet, while 
the rest get transmitted and subsequently relax in the Pt 
layer. The ratio between these two fractions equals the 
ratio between the conductance of the Py-Cu contact and 
the Co-Pt contact, g'^^/g, and is of the order of unity. 
This results in a large magnetization damping as a sig- 
nificant portion of the spin pumping relaxes by spin-orbit 
scattering in Pt. When L is increased, less spins manage 
to diffuse through the entire Cu buffer, and, in the limit 
L ^ Asd, the majority of the spins scatter back into the 
ferromagnet or relax in Cu not feeling the presence of 
the Pt layer at all. In the intermediate regime, the spin 
pumping into the Pt layer has an algebraic fall-off on the 
scale of the elastic mean free path and exponential one 
on the scale of the spin-diffusion length. 

It is important to emphasize that the strong depen- 
dence of damping on the Cu layer thickness L in the Py- 
Cu-Pt configuration gives evidence of the spin accumu- 
lation in the normal-metal system. This spin accumula- 
tion, in turn, indicates that an excited ferromagnet (as in 
the FMR experiment discussed here) transfers spins into 
adjacent nonmagnetic layers, confirming our claimD Fuc^ 
thermore, this supports our concept of the spin battery.E3 

Before ending this section, it is illuminating to make a 
small digression and further study Eq. ( |2^ ) in the limit of 
vanishing spin- flip processes in the buffer layer Nl. Re- 
calling our definitions for Asd and Sgd [Eqs. (0) and (12)] 
and taking limit Tsf ^ oo, we find that Eq. (p6|) reduces 
to Eq. ([2^), only with g^^ replaced by gl^ [similarly to 
Eq. (§)F 



— jTj- — — j-7 + i?jvi H , 



n 



(28) 



V. CONCLUSIONS AND DISCUSSIONS 

Ferromagnets emit a spin current into adjacent nor- 
mal metals when the magnetization direction changes 
with time. We recently proposed a novel mechanism 
for this spin transfer based on the picture of adiabatic 
spin pumping.u It was shown that our theory explains 
the increased magnetization damping in ferromagnets in 
contact witiujiQEmal metals in measurements of FMR 

iinewidths.yy'yEi 

Whereas the spin pumping affects the magnetization 
dynamics, it also creates a nonequilibrium magnetiza- 
tion in adjacent nonmagnetic films. In this paper we 
first calculate this spin accumulation for F-N metallic 
multilayers and find that it induces a spin backflow into 
the ferromagnetic layer that reduces the spin pumping. 
This spin-accumulation-driven current is significant for 
light metals or metals with only s electrons in the conduc- 
tion band, which have a small spin-flip to spin-conserving 
scattering ratio. 

The picture of ferromagnetic spin pumping and subse- 
quent spin diffusion in the adjacent normal- metal layers 
is also applied to the F-NlpN2 configuration in order to 
analyze recent experimentstj on magnetization damping 
in Py-Cu-Pt trilayers. We showed that our theory quan- 
titatively explains the experimental findiugs^ Our analy- 
sis of the experiments by Mizukami et alB\l3 shows that 
FMR of ultrathin ferromagnetic films in contact with sin- 
gle or composite normal-metal buffers is a powerful tool 
to investigate interfacial transport properties of magnetic 
multilayers as well as the spin relaxation parameters of 
the normal-metal layers. 
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where R^i is the resistance of the Nl layer. The right- 
hand side of Eq. ( p8|) is simply the inverse mixing con- 
ductance of the Nl buffer in seriesjvith its two interfaces 
(one with F and one with N2)^ in particular, when 
layer Nl is thick enough, the total mixing conductance 
gig is just the conductance ot-Jihf diffusive normal-metal 
spacer separating F and Af^.Elrc3 The spin pumping into 
layer Nl with the subsequent diffusion and then spin ab- 
sorption by the ideal spin sink N2 (as discussed in this 
section) can thus be viewed as the spin pumping across an 
effective combined scatterer separating the ferromagnet 
(F) from the perfect spin sink (A^^) [as done in obtaining 
Eq. (E3)]. This shows consistency of our approach. 



APPENDIX: ADIABATIC SPIN PUMPING 

Here we present a detailed discussion of spin pumping 
into normal-metal layers by a precessing magnetization 
direction m of an adjacent ferromagnet. A schematic of 
the model is displayed in Fig. ^. The ferromagnetic layer 
F is a a spin-dependent scatterer that governs electron 
transport between two [left (L) and right (i?)] normal- 
metal reservoirs. 

The 2x2 operator /; for the charge and spin current 
in the ^th lead {I = L,R) can be expressed in terms of 
operators aam,i{E) [bam,i{E)] that annihilate a spin-a 
electron with energy E leaving [entering] the Zth lead 
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FIG. 6: Ferromagnetic film (F) sandwiclied between two 
normal-metal layers (N). The latter are taken to be reservoirs 
in common thermal equilibrium. The reflection and transmis- 
sion amplitudes r and t' shown here govern the spin current 
pumped into the right lead. 



^i(E~E')t/h 



through the mth channel: 
if(t)='-Y^ jdEdE' 



(29) 



When the scattering matrix s^'^ {t) of the ferromag- 
netic layer varies slowly on the time scales of electronic re- 
laxation in the system, an adiabatic approximation may 
be used. The annihilation operators for particles entering 
the reservoirs are then related to the operators of the out- 
going states by the instantaneous value of the scattering 
matrix: bam.ii^) = s"nn,wi'^)°'l3n,i'{E)- In terms of aam.i 

only, we can evaluate the expectation value (^i"^{t) 

of the current operator using {al^jn i{E)apn,i'{E')) — 
fiiE)Sai3SmnSwS{E - £"), whcrc fi{E) is the (isotropic) 
distribution function in the Zth reservoir. When the scat- 
tering matrix depends on a single time-dependent param- 
eter X{t), then the Fourier transform of the current ex- 
pectation value Ii{u!) = J dte^'^*Ii{t) can be written as 



Ii{uj) = gx,i{uj)X{uj) 



(30) 



in tftmis of a frequency to- and X-dependent parameter 
9x..iM 



47r ^ 



E 



v 
ds. 



dE 



dfi'jE) 
dE 



dX 



(31) 



Equation ( pO[ ) is the first-order (in frequeaGjy) correc- 
tion to the dc Landauer-Biittiker formula.E2l At equi- 
librium fR{E) = fL{E), Eq. ^ is the lowest-order 



nonvanishing contribution to the current. Furthermore, 
at sufficiently low temperatures, we can approximate 
-~dfi{E)/dE by a 5-function centered at Fermi energy. 
The expectation value of the 2x2 particle-number op- 
erator Qi{lo) [defined by Ii{t) = dQi{t)/dt in time or by 
Ii{lu) = —iLuQi(uj) in frequency domain] for the Ith reser- 
voir is then given by 



\ mnl' 



H. 



\X{lo), (32) 



where the scattering matrices are evaluated at the Fermi 
energy. Because the prefactor on the right-hand side of 
Eq. ( p2[ ) does not depend on frequency lo, the equation is 
also valid in time domain. The change in particle number 
5Qi{t) is proportional to the modulation 5X{t) of param- 
eter X and the 2x2 matrix current (directed into the 
normal-metal leads) reads 



dhi dX{t) 



dX dt ' 
where the "matrix emissivity" into lead I is 
1 ds 



dhi 
dX 



— Y 

mnl' 



mnAl' , TT 



dX 



(33) 



(34) 



If the spin-flip scattering in the ferromagnetic layer is 
disregarded, the scattering matrix s can be written in 
terms of the spin- up and spin-down scattering coefficients 



and 



using the projection matrices 
= (i - ^ • m) /2:EI 



1 



/2 



(35) 



The spin current pumped by the magnetization preces- 
sion is obtained by identifying X{t) — ip{t), where ip is 
the azimuthal angle of the magnetization direction in the 
plane perpendicular to the precession axis. For simplic- 
ity, we assume that the magnetization rotates around the 
y axis: m = (siniyj, 0, cos (^f). Using Eq. (^, it is then 
easy to calculate the emissivity [Eq. (pj)] for this process: 



Qti 1 

— = [ArCTy + Ai(aj; COS ip- sincp)] , (36) 

where Ar{Ai) — Re{lm)[g^ ^ —t^^], as explained in Sec. ||. 
Expanding the 2x2 current into isotropic and traceless 
components. 



-a- Ic 



(37) 



we identify the charge current Ic and spin current Ig. 
Comparing Eqs. (H), @), and (|3|), we find that the 
charge current vanishes, Ic = 0, and the spin current 

h = {Ai cos ip^Ar,- At sin (f) (38) 

47r at 
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can be rewritten as Eq. (|l|). Because the spin current 
transforms as a vector, it is straightforward to show that 
Eq. (]^) is also vahd in the case of the general motion of 
the magnetization direction. 

Even though the mathematics of our scattering ap- 
proach to adiabatic spin pumping is entirely analogous 
to the charge-pumping theory developed in Ref. ^, there 
are some striking differences in the physics. In the case 
of a spin-independent scatterer as in Ref. ^, the average 
charge-pumping current has the same direction in the 
two leads, by charge conservation: the charge entering 
the scattering region through either lead must leave it 
within a period of the external-gate variations. Whereas 



the particle number of the two reservoirs must (on aver- 
age) be conserved also here, the total conduction-electron 
spin angular momentum is not conserved. In fact, as 
we explained in Ref. ^ for a symmetric system shown in 
Fig. ^ a processing ferromagnet loses angular momen- 
tum by polarizing adjacent nonmagnetic conductors. In 
this respect, the phenomenon looks more similar to a 
spin "well" or "fountain" : An excited ferromagnet ejects 
spins in all directions into adjacent conductors by los- 
ing its own angular momentum, rather than transfers 
("pumps") spins from one lead to the other. The an- 
gular momentum has to be provided, of course, by the 
applied magnetic field. 
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